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Abstract. Spin Hall effect converts charge current to spin current, which can exert 
spin-torque to switch the magnetization of a nanomagnet. Recently, it is shown that 
the ratio of spin current to charge current using spin Hall effect can be made more 
than unity by using the areal geometry judiciously, unlike the case of conventional 
spin-transfer-torque switching of nanomagnets. This can enable energy-efficient means 
to write a bit of information in nanomagnets. Here, we study the energy dissipation 
in such spin Hall devices. By solving stochastic Landau-Lifshitz-Gilbert equation of 
magnetization dynamics in the presence of room temperature thermal fluctuations, we 
show a methodology to simultaneously reduce switching delay, its variance and energy 
dissipation, while lateral dimensions of the spin Hall devices are scaled down. 
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1. Introduction 

Spintronics is a promising field of research that can possibly replace the traditional 
charge-based electronics HI El Spin-transfer-torque (STT) is a current-induced 
magnetization switching (CIMS) mechanism in which magnetization of a nanomagnet 
can be switched between its two stable states [a El El E| or it may lead to 
oscillatory motion of magnetization [8], however, the energy dissipation in such 
switching mechanism is too high for practical application purposes compared to 
the traditional charge-based electronics [9]. There are other mechanisms that are 
coming along for energy-efficient switching of a bit of information e.g., electric field- 
induced magnetization switching in multiferroic heterostructures [101 IHl 1121 ES ITT] , 
perpendicular anisotropy [151 ESI EZ] , coupled polarization-magnetization switching in 
single-phase multiferroic materials HHl ESI Efll EH- However, one mechanism that has 
gotten attention recently is to utilize the spin Hall effect, which was hrst recognized 
by D’yakonov and Perel’ [22| and following which there have been both theoretical 
studies [23l 1211 [25l 126] and experimental investigations in semiconductors [271 EH ESI [30] . 
Utilizing spin Hall effect to generate a sufficient spin current for technological purposes 
of exciting magnetization dynamics was severely limited Ell. however, there have been 
recent resurgence of interests [311 |32l |33l [3l| due to giant spin Hall effect of exerting spin- 
torque, which can be used to switch the magnetization direction of a nanomagnet using 
different spin Hall materials e.g., platinum [35l ESI [371 [38] , tantalum [39], tungsten [lO] . 
CuBi [IT] . 

Fig. El shows a schematic diagram of a spin Hall device. A charge current flows 
laterally in the spin Hall material layer and a perpendicular spin current through the 
cross-section of the structure is generated, which can exert spin-torque on the free layer 
nanomagnet and switch its magnetization. In the conventional STT devices, the charge 
current gets spin polarized through a hxed layer, however, passing high current during 
switching of magnetization occasionally damages the thin insulator layer [39]. Hence, 
the giant spin Hall devices offer a great advantage over the traditional STT devices. 
Note that the read current, which is small, still needs to be passed perpendicular to 
the cross-section of the spin Hall devices and hence the fixed layer in a magnetic tunnel 
junction (MTJ) structure is required to detect the magnetization states of the free layer 
nanomagnet [12]. By reversing the direction of the charge current, the magnetization 
of the free layer can be switched in either direction. Since the charge current and the 
spin current flow though different cross-sections, it is shown that the spin-to-charge 
current ratio can be greater than one, unlike the case of traditional STT devices [39]. 
Both the Oersted field generated by the current and Rashba held are neglected since 
experimentally it has been verihed that they tend to oppose the switching and of small 
magnitude [39] . 

Here, we perform an analysis on the energy dissipation in the aforesaid spin Hall 
devices. Since the charge current hows through the spin Hall material, the Ohmic loss 
occurs in the spin Hall material itself unlike the case of traditional STT devices, where 
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Figure 1. Simplified schematic diagram of a giant spin Hall device and axis 
assignment. The charge current (of current density Jc) flows in the y-direction through 
the spin Hall material layer accumulating spins of opposite polarities on its opposite 
surfaces, which generates a spin current (of current density Jg) in the ^-direction 
with spin polarization in the z-direction. For the direction of charge current and spin 
current shown in the diagram, the direction of spin polarization at the free layer-spin 
Hall material layer interface will be in the -|-z-direction for materials with positive spin 
Hall angle (e.g., Pt), while it will be in the —z-direction for materials with negative 
spin Hall angle (e.g., CuBi). Since the two mutually anti-parallel stable magnetization 
orientations of the free layer along the z-axis encode the logic bits 0 and 1, the spin- 
torque acting in the free layer needs to switch the magnetization in either direction and 
it can be performed by changing the direction of the charge current. Note that unlike 
the case of traditional spin-transfer-torque devices, the charge current is not spin- 
polarized by the fixed layer, however, the fixed layer and the insulator (e.g., MgO) are 
required to read the magnetization state of the free layer by passing a small charge 
current in the a;-direction via tunneling magnetoresistance (TMR) mechanism. In 
standard spherical coordinate system, 6 is the polar angle and (/) is the azimuthal 
angle. 


the energy dissipation incnrs in the MTJ stack. The charge cnrrent flows laterally 
throngh a higher resistance compared to the case if flown perpendicnlar to the cross- 
section (of the nanomagnet), however, the MTJ stack (nanomagnets separated by an 
insnlator) has a higher resistance than that of the cnrrent flow path in the spin Hall 
material layer, particularly while using a low-resistivity material for the spin Hall layer. 
Also, the required charge current is less due to higher spin-to-charge conversion factor 
in giant spin Hall devices. We show that the energy dissipation is inversely proportional 
to the lateral area of the spin Hall devices, but scaling down the lateral dimensions also 
makes the required spin current to switch the magnetization smaller. We solve stochastic 
Landau-Lifshitz-Gilbert (LLG) equation P31ITT|, IT5] of magnetization dynamics to depict 
a strategy whereby the energy dissipation can be reduced simultaneously with lowering 
of the switching delay and its variance in the presence of room temperature thermal 
fluctuations, while scaling down the lateral dimensions. 
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2. Model 

We consider the free layer nanomagnet in the shape of an elliptical cylinder with its 
elliptical cross section lying on the y-z plane; the major axis is along the 2 ;-direction 
and the minor axis is along the ^/-direction (see Fig. 1). The dimensions along the z-, 
y- , and x-direction are a, b, and / (a > 6 S> /), respectively. So the lateral area of 
the spin Hall device is H = (7r/4)a6 and the nanomagnet’s volnme is H = (7r/4)a6/. In 
standard spherical coordinate system, 6 is the polar angle and the 0 is the azimuthal 
angle. Due to shape anisotropy, the two degenerate magnetization states along the 2 ;- 
direction {6 = 0° and 180°) can store a binary bit of information. The y-axis is the 
in-plane hard axis and x-axis is the out-of-plane hard axis. We can write the shape 
anisotropy energy of the nanomagnet as follows: 

Eshape{0, (j)) = ^M[Hk + Hd cosV] ^ sin'^e, (1) 

where M = fioMg, po is permeability of free space, is the saturation magnetization, 
Hk = {Ndy—Ndz) Mg is the Stoner-Wohlfarth switching held [16], Hd = (Ndx—Ndy) Mg is 
the out-of-plane demagnetization held [IB], and Ndm is the m*^ (m = x, y, z) component 
of the demagnetization factor {Ndx 3> Ndy > Ndz)- The shape anisotropic energy 
barrier between the two stable magnetization states can be expressed from Equation (|Ti) 
[putting 0 = ±90°] as 

Ebarrier = ^M Hkfl. ( 2 ) 

The magnetization M of the nanomagnet has a constant magnitude but a variable 
direction, and hence we can represent it by a vector of unit norm n„i = M/|M| = br, 
where br is the unit vector in the radial direction in spherical coordinate system 
represented by (r,6*,0); the other two unit vectors in the spherical coordinate system 
are Bq and b^ for 9 and 0 rotations, respectively. 

The ehective held and torque acting on the magnetization due to gradient 
of potential landscape can be expressed as Heff(6*,0) = —VEghape{G) 4 ') = 

-{dEghape/d 6 )e 0 - (l/sine) {dEghape/d(j)) and Te(6',0) = x Heff (0,0), 

respectively. 

0) Bgfiape(^4^') siTl(2,9'^ B(j) Bgiiape,(f>{_4^^ sin9 Gg, (3) 

where 

Bghapeifp) = (1/2) M [Hk + Hd cosV] (4) 

BghapeA^) = (1/2) MifrfDsm(2(/). (5) 

Passage of a spin current Ig in the x-direction with spin polarization along the 
2 ;-direction generates a spin-transfer-torque that is given by [4] 

TsTTit) = -sn^{t) X (iTs X nm(f)) = ssinOeg, (6) 

where s = {h/ 2 e)Ig is the spin angular momentum deposition per unit time, the 
unit vector rig = ±bz is in the direction of spin polarization since we want to 
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rotate magnetization from —z-axis to +z-axis, and we have used the identity = 
cos6 hr — sin6 Bq. No asymmetry in spin-transfer-torque switching or held-like torque is 
considered here [IHl HHl [56]. 

The random thermal fluctuations are incorporated via a random magnetic held 
h{t) = hx{t)ey^ + hy{t)ey + hzit)ez, where hi{t) {i = x, y, z) are the three components of 
the random thermal held in Cartesian coordinates. We assume the same properties of 
the random held h(f) as described in Ref. |15]. The random thermal held can be written 
as [15] 


hi(t) 


2akT 

\-f\MnAt 


G(o,i)(t) {iex,y,z) 


(7) 


where a is the phenomenological dimensionless Gilbert damping parameter, k is the 
Boltzmann constant, T is temperature, 7 is the gyromagnetic ratio for electrons and its 
magnitude is equal to 2.21 x 10 ® (rad.m).(A.s)“^, At is the simulation time-step used, 
and the quantity G(o,i)(t) is a Gaussian distribution with zero mean and unit variance. 

The thermal held and the corresponding torque acting on the magnetization 
can be written as Hth( 6',0, t) = Pe{9, bq + P^{9, (j),t) b^ and Tth( 6',0, t) = 
rim X 0,^)5 respectively, where 


Fe(6*, 0, t) = M[hx{t) cos9 coscj) + hy{t) cos9sin(j) — hz{t) sin9] G, (8) 

P^{9, (j), t) = M[hy{t) cos(j) — hx(t) sm0] G. (9) 

The magnetization dynamics under the action of the torques TE(t), Tstt(^), and 
Tth(^) is described by the stochastic Landau-Lifshitz-Gilbert (LLG) equation jlSUHl 
05] as follows. 

From the above equation, we get the following coupled equations of magnetization 
dynamics for 9 and 0: 

^ ^ ^ - s) sin9 

- aBshapei<P)sin{29) + {aP0{9, 0, t) + P^{9, 0, t))], (11) 


^ ^ ^ -s) + 2B{(j))cos9 

— [sin9]~^ {Pe{9, 0, t) — aP^{9, 0 , f))] {sin9 7 ^ 0 ). ( 12 ) 

We solve the above two coupled equations numerically to track the trajectory of 
magnetization over time. 

The energy dissipated in the nanomagnet due to Gilbert damping can be expressed 
as Edamp = Jq Pdampit)dt, where r is the switching delay and Pdampif) is the power 
dissipated at time t given by 

Pdampit) = |TE( 6 '(t), 0(t))T stt(^)P- 


(13) 
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Thermal field with zero mean does not canse any net energy dissipation bnt it canses 
variability in the energy dissipation by scnttling the trajectory of magnetization. 

If the magnetization sitnates exactly along the easy axis, i.e., sin^ = 0 (0 = 0° 
or 6 = 180°), the torqne acting on the magnetization given by Eqnations ([3]) and (]6|) 
becomes zero. That is why only thermal finctnations can deflect the magnetization 
vector exactly from the easy axis. Magnetization finctnates aronnd an easy axis dne to 
thermal agitations and hence we get a distribntion of the initial angles [9 inuiah^t)initial)- 
We consider this initial distribntion when magnetization starts switching from 9 ~ 180°. 
We perform a moderately large nnmber (10,000) of simnlations in the presence of thermal 
finctnations and when the final valne 9final becomes < 4.5° (note that the mean valne 
of 9iniuai is ~4.5°), the switching is deemed to have completed. Then, the mean and 
standard deviation of switching delay distribntion {xmean and Tstd, respectively), and the 
mean of energy dissipation Edamp,mean are extracted from the simnlations. (See Fig. [2] 
and Table [T] later.) 

The energy dissipation dne to Gilbert damping is of the order of the energy barrier 
height, however, the major part of the energy dissipation occnrs dne to Ohmic loss in 
the spin Hall material since the charge cnrrent Ic flows throngh it in the ?/-direction (see 
Fig. [1]). The spin-cnrrent Ig flows in the x-direction, and the ratio of the spin-to-charge 
cnrrent can be expressed as 


where Jg (Jc) is the spin (charge) cnrrent density, Ag (Ac) is the area throngh which 
spin (charge) current flows, Qsh is the spin Hall angle of the spin Hall material used, 
and tsH is the thickness of the spin Hall material layer. The geometric factor b/tsn can 
be selected much greater than one and hence, having Qsh b > tsn can make Ig > Ic- 
The power dissipation in the spin Hall material can be expressed as 



(15) 


where p is the resistivity of the spin Hall material and R is the resistance of the spin Hall 
material layer. The energy dissipation in the spin Hall material layer can be expressed 
as 



(16) 


where Tp is the time-period until which the charge current is kept active. 


3. Results 


For the free layer nanomagnet, we consider the widely-used material CoFeB, which 
has low saturation magnetization Mg = 8 x 10^ A/m [50] and a low Gilbert damping 
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Figure 2. Switching delay distributions from 10,000 simulations by solving stochastic 
Landau-Lifshitz-Gilbert (LLG) equation of magnetization dynamics in the presence 
of room temperature (300 K) thermal fluctuations, (a) a=150 nm, 6=100 nm, 1=2 
nm, and Ic = 133.3 /rA. The mean of switching delay distribution and its standard 
deviation are 1.03 ns and 0.21 ns, respectively, (b) a=100 nm, 6=30 nm, 1=2 nm, and 
Ic = 198.8 /iA. The mean of switching delay distribution and its standard deviation 
are 0.55 ns and 0.10 ns, respectively, (c) a=100 nm, 6=30 nm, 1=2 nm, and Ic = 133.3 
/rA. The mean of switching delay distribution and its standard deviation are 0.78 ns 
and 0.15 ns, respectively. The distributions look like log-normal distributions and can 
be fitted accordingly. The means and the standard deviations are directly calculated 
from the numerically computed distributions. 


parameter a = 0.01. Note that the damping parameter can be modihed depending on the 
adjacent spin Hall material layer due to spin pumping [51]. We choose CuBi as the spin 
Hall material since it has the lowest p/Q'sh factor (p=10 fiQ-cm, 05'jy=-O.24 measured 
at 10 K |4T]) among the other materials utilized in current experiments with an eye to 
reduce the energy dissipation [see Equation (ITOl) ]. We can also possibly utilize platinum 
as the spin Hall material layer that has p ~ 15 pVt-cm [36], 05 jy=O.O 7 [38] measured 
at room temperature, however, it will incur an energy dissipation of one order more 
than that of CuBi. Note that platinum increases magnetization damping (i.e., increases 
the minimum switching current) of CoFeB considerably and the damping parameter is 
about three times compared to when tantalum is utilized [39]. 

Note that CoFeB has a resistivity of 100 pQ-cm [52], which is 10 times more than 
that of CuBi spin Hall material layer. Hence the current shunting effect [53] IM] through 
the CoFeB layer is ignored. Therefore, the purpose of utilizing a low-resistive spin 
Hall material layer is also to tackle the issue of current shunting effect, apart from 
reducing the energy dissipation. We do need to consider the current shunting effect 
while utilizing tantalum or tungsten as the spin Hall material layer since they have 
the resistivity of the same order as of CoFeB [391 HO]. We consider here in-plane 
switching of magnetization, however, perpendicular switching of magnetization has been 
demonstrated [TH] [55] EE] 139] and can be considered too, particularly to scale down the 
lateral area of the devices further. 

We use both the thicknesses tsH and I as 2 nm. We vary the lateral dimensions a 
and b of the nanomagnet for different cases considered and we choose the dimensions 
of the nanomagnet so that it has a single ferromagnetic domain [53 HZ]- always 


















Table 1. Different metrics for three different cases considered. Cases (a) and (b) 
correspond to the same dissipated power Pd= 377 nW, while the cases (a) and (c) 
correspond to the same charge current Ic = 133.3 /iA. The energy dissipation is reduced 
by about 5 times for the case (c) compared to the case (a), while the lateral area is 
scaled down by 5 times. 


Case 

Nanomagnet size 
(nm3) 

Is 

(mA) 

Ic 

(mA) 

R 

m 

'^mean 

(ns) 

'^std 

(ns) 

Ej damp,mean 

(aJ) 

Tp 

(ns) 

Pd 

(nW) 

Ed 

(aJ) 

(a) 

150 X 100 X 2 

1.6000 

0.1333 

21.22 

1.03 

0.21 

0.0883 

2.5 

377 

943 

(b) 

100 X 30 X 2 

0.7155 

0.1988 

9.55 

0.55 

0.10 

0.0395 

1.2 

377 

453 

(c) 

100 X 30 X 2 

0.4800 

0.1333 

9.55 

0.78 

0.15 

0.0485 

2.0 

170 

204 


keep the energy barrier Ebarrier = 0.8 eV or 32 kT at room temperature (T=300 
K). Following Boltzmann distribution, this means that the error probability due to 
spontaneous reversal of magnetization is or 10“^^, which is low 

enough for application purposes. We can choose a higher barrier height leading to even 
lower error probability depending on the application requirements. We will now go 
through the following three cases. 

Case (a): We consider a nanomagnet with dimensions a = 150 nm, b = 100 nm, 
and I = 2 nm. We solve stochastic LLG for 10,000 times and hnd that it requires about 
Is = 1.6 mA to switch always successfully, while the switching current is kept active 
for 2.5 ns. The corresponding charge current is A = 133.3 pA. The demagnetization 
factor {Nd:,, Ndy, Nd,) = (0.9468, 0.0339, 0.019) |17], Hk = 0.0148 T, and Hd = 0.9177 
T. The switching delay distribution is plotted in the Fig. [2](a). The mean of switching 
delay and its standard deviation are 1.03 ns and 0.21 ns, respectively. The major energy 
dissipation is due to Ohmic loss in the spin Hall material and it turns out to be 943 aJ 
[see Table [I]. 

Case (b): We now reduce the lateral dimensions of the nanomagnet and choose the 
lateral dimensions to be a = 100 nm, 6 = 30 nm with an eye to increase the device 
density on a chip. The thickness I is kept same as for the case (a). Hence, both the 
lateral area and the volume decrease by a factor of 5. Note that the energy barrier 
between the two stable states is kept constant [see Equation (|2])] at this reduced volume 
by choosing the elliptical cross-section more anisotropic. The demagnetization factor 
(A'dx, Ndy, Ndz) = (0.8829, 0.0953, 0.0218) Hk = 0.0739 T, and Hd = 0.7917 T. 
The anisotropic held Hk increases due to the modihcation of demagnetization factors at 
the chosen dimensions m- However, at a reduced volume, the magnetization becomes 
more prone to thermal huctuations [see Equation ([7])]. But the STT switching current 
also mitigates the detrimental effects of thermal agitations. Since in this case the current 
needs to switch a nanomagnet of a smaller volume [see Equations flTTll and ffT2ll ]. we 
show that it is possible to adjust the switching current such that in overall, we can 
achieve a better performance metrics for both switching delay and energy dissipation at 
this reduced volume. 
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According to the Equation ffTBD . we scale down A by a factor of \/5 to keep the 
power dissipation constant. The corresponding distribution of switching delay is shown 
in the Fig. [2]^b). Note that both the mean and standard deviation of switching delay 
have got reduced by half compared to the case (a). The switching current is kept on 
until 1.2 ns and hence the energy dissipation is reduced by more than half compared 
to the case (a). Note that even if the spin current A is reduced, the charge current 
required in this case has got increased compared to the case (a) due to the decrease in 
minor axis b of the elliptical cross-section of the nanomagnet [see Equation (ITTll ]. Next 
we consider another case where the spin current A is reduced further to keep the charge 
current same compared to the case (a). 

Case (c): We choose the same dimensions of the nanomagnet as of the case (b) 
and the same charge current as for the case (a) [see Table [1]. The demagnetization 
factor, TTfc, and are same as of the previous case. The corresponding switching 
delay distribution is plotted in the Fig. |2][c). Note that both the mean and standard 
deviation of switching delay distribution have got increased compared to the case (b) 
due to decrease in A, however, the metrics are still much better than that of the case (a) 
[see Table [T] . The energy dissipation Ed has got reduced by more than half, compared 
to the case (b), to 204 aJ, which is due to the decrease in spin current Ig. 

It should be noted that it is not only the mean of switching delay distribution, but 
also the standard deviation in switching delay that plays an important role in setting 
the clock period of magnetization switching for application purposes. A higher standard 
deviation would lead to setting a higher clock period. In this paper, it is shown that both 
the mean and the standard deviation of switching delay can be reduced simultaneously, 
while the lateral area of the giant spin Hall devices is scaled down. 

4. Discussions 

From Equation (ITTjl . it should be noted that the spin diffusion length \sh for spin 
Hall material is not considered in the expression. For a thick spin Hall material layer 
(tsH ^sh), if does not need to consider Xsh- However, we should choose the thickness 
of the spin Hall material layer tsH small since the charge current (for a required spin 
current) decreases with the decrease of tsH- If is possible to add a spin-sink layer at the 
bottom of a thin spin Hall material layer (see Fig. [1]) to avoid any backflow of spins. For 
example, see Ref. [58] for some experimental results, however, research on such front is 
quite emerging. For CuBi utilized as a spin Hall material layer, such experimental data 
is not available, however, this concept of adding a spin-sink layer is quite general. Also it 
should be noted that there is controversy on the experimentally measured spin diffusion 
length for spin Hall materials [5H]- Since for all the three cases in Tabled] the results 
correspond to the same thickness of the spin Hall material layer, the comparative nature 
of the analysis is not quite affected. The analysis presented here depicts the necessity of 
decreasing the thickness of the spin Hall material layer to decrease the charge current 
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and consequently to reduce the energy dissipation. 

It is imperative to compare the energy savings utilizing the giant spin Hall effect 
compared to the traditional way of exerting spin-torque. The energy savings accrue from 
the decrease in charge current due to the geometric factor and decrease in resistance of 
the charge flow path compared to the MTJ stack. The reduction of energy dissipation 
is as high as 3-4 orders of magnitude compared to the conventional spin-transfer-torque 
switching mechanism [59l [60] and domain-wall racetrack memory [611 [62] . 

Although the energy dissipation in these giant spin Hall devices has got reduced to 
the order of 0.1 fJ, the energy dissipation can be further reduced by having a spin Hall 
material that has even lower factor [see Equation flTOll ]. i.e., having a material 

with a lower resistivity and a higher spin Hall angle. The target would be the reduction 
of energy dissipation by a factor of 2 more to be competitive with other emerging 
technologies [10] . Switching delay and area of a device also need to be competitive with 
the traditional transistor based technology [63]. The non-volatility of the nanomagnets 
in these giant spin Hall devices can be utilized to devise a possibly better architecture 
in terms of performance metrics for application purposes. 

5. Conclusions 

In conclusion, we have analyzed the energy dissipation in recently proposed spin Hall 
devices exploiting giant spin Hall effect. After formulating the energy dissipation, we 
solved stochastic Landau-Lifshitz-Gilbert equation of magnetization dynamics in the 
presence of room-temperature thermal fluctuations to present a methodology in which 
both the energy dissipation and switching delay (and its variance) can be reduced 
simultaneously, while the lateral dimensions of the spin Hall devices are scaled down. 
The energy dissipation turns out to be of several orders of magnitude less than that of 
the traditional spin-transfer-torque devices. This held is still emerging and with suitable 
spin Hall materials, the energy dissipation can be reduced further. This opens up an 
energy-efficient avenue to control a bit of information in nanomagnetic memory and 
logic systems for our future information processing paradigm. 
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